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Abstract. Quantum error correction provides a fertile context for exploring the 
interplay of feedback control, microscopic physics and noncommutative probability. 
In this paper we deepen our understanding of this nexus through high-level analysis of 
a class of quantum memory models that we have previously proposed, which implement 
continuous-time versions of well-known stabilizer codes in autonomous nanophotonic 
circuits that require no external clocking or control. We show that the presence 
of the gauge subsystem in the nine-qubit Bacon-Shor code allows for an optimized 
layout of the corresponding nanophotonic circuit that substantially ameliorates the 
effects of optical propagation losses, argue that code separability allows for simplified 
restoration feedback protocols, and propose a modified fidelity metric for quantifying 
the performance of realistic quantum memories. Our treatment of these topics exploits 
the homogeneous modeling framework of autonomous nanophotonic circuits, but the 
key ideas translate to the traditional setting of discrete time, measurement-based 
quantum error correction. 
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Quantum feedback control provides a systems engineering perspective on the analysis 
and design of quantum memories, complementing alternative ideas extending from 
theoretical physics. Whereas the latter approach has emphasized the connections of 
quantum error correction and quantum computation to many body physics (TJ |2], the 
quantum control community has viewed decoherence suppression as a problem that 
should ideally be formulated as a non-commutative generalization of classical stochastic 
and hybrid control theory. In a sense, we can view this research program as a 1950's/60's- 
era agenda translated forward a half century. That is, given the novel resources of 
photonics, quantum electronics, and spintronics, we can ask ourselves what the basic 
constituent components might be of a post-classical information processing machine, and 
how the exchange of signals among such components could be used to realize desired 
functionality [31 HIE]. 

Within this setting we have recently investigated nanophotonic circuit models for 
quantum memories that autonomously implement well-known stabilizer codes. In these 
models, each physical qubit is strongly coupled to an optical or microwave resonator, 
and the resonators are coupled by waveguides to form a coherent feedback network 
(photonic circuit). When the circuit is powered by appropriate stationary laser inputs, 
the resulting continuous-time dynamics effect error detection and correction without 
any additional external clocking or control signals. In such models, the classical signal 
processing apparatus assumed in standard treatments of quantum error correction 
is replaced by a small number of controller qubits, making it possible to derive a 
master equation for the closed-loop behavior that can be modified straightforwardly 
to incorporate a wide range of realistic decoherence mechanisms [6l [TJ, [8] . 

The operational principles of these quantum memory models, as well as the methods 
used to derive their equations of motion, exemplify a quantum-optical generalization 
of conventional electric circuit theory in which guided electromagnetic fields play the 
role of signals and qubit-resonators serve as input /output devices that process them 
coherently. We have previously shown how the the Gough-James quantum network 
algebra [HJ QUI [EE] can be utilized together with a recent limit theorem for quantum 
stochastic differential equations (QSDEs) [12J to facilitate the derivation of an intuitive 
master equation for a given quantum memory model from an explicit construction 
of the underlying nanophotonic circuit, in a manner inspired by schematic capture 
methods of contemporary electrical engineering [8j [13]. Here, in order to focus on 
higher architectural principles, we will skip over such details and jump directly to 
considering a class of master equations that arises from the general approach of 
embedding continuous-time relaxations of stabilizer quantum error correcting codes in 
the autonomous dynamics of a coherent feedback network. 

We model an autonomous quantum memory using a Markovian master 
equation [TJ], 



where if is a Hamiltonian for the internal dynamics of the memory and controller 




(1) 
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qubits, and the Lindblad operators describe couplings to reservoir modes that mediate 
decoherence (including memory errors) and entropy removal (via spontaneous emission- 
type processes). At this level of description, the dynamics of the electromagnetic 
field modes within the waveguides and resonators of the photonic circuit have been 
adiabatically eliminated. Details of the derivation of such a master equation for the 
three-qubit code can be found in |6]; here in considering arbitrary stabilizer codes we 
assume that the register qubit model [TH] can be generalized to incorporate atomic level 
structures whose Raman resonance conditions effectively implement an AND operation 
on as many feedback signals as the code requires. We assume a feedback controller 
construction based on nanophotonic relays [H] as in prior work. 

Consider a stabilizer quantum error correcting code [201 EI] that encodes one qubit 
of information in Q register qubits with N stabilizer generators. We refer to the complete 
set of stabilizer generators as S = {M n }^ =1 . Recall that a stabilizer code redundantly 
encodes information in the joint state of the register qubits, that each stabilizer generator 
is a joint observable of the register qubits with eigenvalues ±1, and that a measurement 
of the full set of stabilizer generators values (the error syndrome) suffices to localize 
any correctable error but yields no information on the encoded qubit. To explicitly 
construct the master equation for an autonomous quantum memory based on S we 
use the following procedure. For each stabilizer generator M n include two Lindblad 
operators of the form 

L 2n _! = a (a** (I + M n ) - nf (I - M n )) , (2) 
L 2n =a(o- R »(I-M n ) + Tl^(I + M n )), (3) 

where erf™, H Rn and IT^ n are (qubit-like) raising/lowering operators and projectors onto 
g and h basis states for the n th relay of the feedback controller, and a parameterizes 
the amplitudes of the electromagnetic probe fields, which in turn determines the 
syndrome measurement rate. If the code is separable, in the sense that disjoint sets of 
stabilizer generators mediate X and Z error syndrome extraction, the quantum memory 
Hamiltonian can be written in the form 
Q 

H = Y,^{X n - F Sx [X n \ + Z n ■ T Sz [Z n \ ) . (4) 

n=l 

If the code is not separable then 
Q 

H = J2^( X n-J r s [X n ] +Y n -7 S [Y n ] + Z n -7 S [Z n \ ) , (5) 

71=1 

where the function Fs [E n ] maps a single-qubit error to a projector onto a state of 
the controller relays. The action of J 7 straightforwardly represents the lookup table of 
correctable errors and corresponding syndromes. For example, suppose there are six 
stabilizer generators and that they take the values 1-1-11-11 for a register state with 
a bit-flip error on the fourth qubit. Then JF S [X 4 ] = U Rl Uf 2 Ilf 3 Il R4 Ilf 5 U R6 (the h state 
of relay n is associated with the +1 value of M n , and the g state with value —1). When 



Gauge subsystems, separability, and robustness in autonomous quantum memories 4 



the code is separable, we use Sx and Sz to denote disjoint subsets of S. In both of the 
above Hamiltonian expressions, Q parameterizes the strength of the feedback fields. 

To illustrate our construction we first consider the seven-qubit code [201 El] , which 
is separable. The following lookup table gives the corresponding error syndromes for 
single-qubit X, Z, and Y errors with ■ corresponding to a +1 value for a given stabilizer 
generator and Q to —1: 



Error syndromes for 7-qubit code 
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It is easily seen that the error syndrome for Y n is simply the logical OR of the syndromes 
for X n and Z n , which combined with the algebraic relation Y n oc X n Z n makes it possible 
to design a feedback network that independently detects and corrects X and Z errors — 
when a Y n error occurs it can be diagnosed and treated as the occurrence of both an 
X n and a Z n error. 

Applying the construction described above for the separable case, we arrive at the 
following master equation for a coherent feedback implementation of the seven-qubit 
code: 

h = n^nf nf nf + x 2 nf nf nf + x 3 nf nf nf + 

x 4 nf nf nf + x 5 nf nf nf + x 6 nf nf nf + 
x 7 nf nf nf + z x nf nf nf + z 2 nf nf nf + 
z :l nf ii,';Tif + z 4 nf nf nf + z 5 nf nf nf + 

z 6 nf nf nf + z 7 nf nf nf ) , (6) 
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L x = a (af (7 + X]X 2 X 3 X 4 ) - Ilf (7 - X 1 X 2 X 3 X 4 )) , 

L 2 = a (a Rl (7 - X!X 2 X 3 X 4 ) + nf 1 (/ + X 1 X 2 X 3 X 4 )) , 

L 3 = a (af (7 + X 1 X 2 X 5 X 6 ) - nf (7 - X^XgXg)) , 

L 4 = a (a* 2 (7 - X 1 X 2 X 5 X 6 ) + nf (7 + X X X 2 X 5 X 6 )) , 

L 5 = a (af (7 + X 1 X 3 X 5 X 7 ) - nf (7 - X 1 X 3 X 5 X 7 )) , 

L 6 = a (a* 3 (7 - X 1 X 3 X 5 X 7 ) + nf (7 + X 1 X 3 X 5 X 7 )) , 

L 7 = a (af (7 + Z X Z 2 Z 3 Z 4 ) - nf (7 - Z x Z 2 Z 3 Zt)) , 

L 8 = a (af (7 - ZiZ 2 Z 3 Z 4 ) + nf (7 + Z X Z 2 Z^Z 4 )) , 

L 9 = a (af (7 + ZxZ 2 Z 5 Z 6 ) - nf (7 - Z^s^)) , 

L10 = « (a^(7 - Z^^e) + nf (7 + Z^^e)) , 

L u = a (af (7 + Z 1 Z 3 Z S Z 7 ) - nf (7 - Z 1 Z 3 Z S Z 7 )) , 

£12 = a {o- R °(I - Z x Z 3 Z h Z 7 ) + nf (7 + Z 1 Z 3 Z 5 Z 7 )) . (7) 

Any desired Markovian error model can be incorporated via additional Lindblad terms, 
e.g., for bit-flip errors, 

L 12+n = VfX n , nel...Q, (8) 

or for spontaneous emission-type decoherence, 

L 12+n = Vf(X n -iY n ), nel...Q. (9) 

In either case, the parameter T adjusts the decoherence rate. Noise processes acting on 
the controller degrees of freedom can be included in the analogous fashion. 

If we next consider the five-qubit code, the smallest quantum error correcting code 
capable of protecting a single encoded qubit against arbitrary single-qubit errors [201 EI] , 
the following lookup table gives the error syndromes for X, Z and Y errors acting on 
the register qubits: 
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It appears by inspection that the syndromes of X n , Z n and Y n are not simply related, 
so evidently we must implement a feedback controller with sufficient logic to react 
conditionally to each of the fifteen four-bit syndromes. For the separable seven-qubit 
code we have seen that only two independent sets of seven three-bit syndromes need 
to be interpreted, suggesting that the feedback control sub-circuit in an autonomous 
quantum memory based on the seven-qubit code could be substantially simpler than for 
the five-qubit code, even though it involves two additional stabilizer generators. Given 
that we have in mind a homogeneous implementation paradigm, in which the quantum 
error- correcting controller is constructed from components that are very similar in nature 
to those of the codeword register [6], this comparison suggests a general advantage of 
separable codes in terms of implementation circuit complexity and concomitant physical 
resource requirements. 

Applying the construction described above for a non-separable code, we arrive at 
the following master equation for a coherent feedback implementation of the five-qubit 
code: 

h = n^nf nf nf nf + x 2 nf nf nf nf + x 3 nf nf nf nf + 

x 4 u R 4 nf nf nf + x 5 nf x nf nf nf + 

^nf nf nf nf + z 2 nf nf nf nf + z 3 nf nf nf nf + 

z 4 nf nf nf nf + z B nf nf nf nf 4 + 

vinf nf nf nf + r 2 nf nf nf nf + y 3 nf nf nf nf + 

r 4 nf nf nf nf + r 5 nf nf nf nf ), (10) 

L 1 = a (cxf (7 + Z 2 X 3 X 4 Z 5 ) - nf (7 - Z 2 X 3 X 4 Z 5 )) , 
7 2 = a (a R i(I - Z 2 X 3 X 4 Z 5 ) + nf (/ + Z 2 X 3 X 4 Z b )) , 
L 3 = a (af (7 + ^^3X4X5) - nf (J - ^3X4X5)) , 
L 4 = a (a^(I - Z 1 Z 3 X 4 X 5 ) + nf (/ + ^^3X4X5)) , 
7 5 = a (af (7 + X!Z 2 Z 4 X 5 ) - nf (7 - X x Z 2 Z 4 X h )) , 
L 6 = a (a R 3(I - X^ZtXs) + nf (7 + X x Z 2 Z 4 X h )) , 
7 7 = a (af (7 + X^^) - nf (7 - X^^)) , 
L 8 = a (a R *(I - X^^) + nf (7 + X x X 2 Z 3 Z h )) . (11) 

Again, decoherence processes acting on the register qubits and/or the controller degrees 
of freedom can be incorporated using additional Lindblad terms. 

The modeling approach that we have described also admits a straightforward 
extension to incorporate the effects of optical propagation losses in the waveguides 
that connect components within the nanophotonic circuit [8], which lead to important 
considerations of optimal circuit layout that will be discussed below. As an illustrative 
example we first consider the probe subnetwork that extracts the error syndrome for 
the first stabilizer generator in the five-qubit code, Z 2 X 3 X 4 Z 5 . This measurement is 
implemented by sequentially interrogating the second, third, fourth and fifth register 
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qubit-resonator components with a coherent laser field. In an idealized model with 
no optical propagation losses for the probe field, the net effect of these couplings is 
represented completely by the pair of coupling terms L\ and L2 given in the above 
section. In a more realistic model that accounts for optical waveguide losses, however, 
information leaks out into the environment as the probe field propagates between 
components in the photonic circuit. It can be shown j8] that the resulting decoherence 
processes are described by additional coupling terms which amount to errors, some of 
which are correlated errors of multiple qubits that the five-bit code is not designed to 
correct: 

L n = a9Z 5 , L n+ i = a9X 4 Z 5 , 

L n +2 = oi9X 3 XiZ 5 , L n+3 = a9Z 2 X 3 X 4 Z 5 . (12) 

Here 9 parameterizes the probe field loss per waveguide segment. Similar sets of 
additional errors would arise from losses along the probe field paths associated with 
each of the other three stabilizer generators. 



While our analysis leading to Eqs. (12) has been grounded in the specialized 



modeling framework of autonomous nanophotonic circuits, our findings generally parallel 
known results from discrete-time, measurement-based implementation scenarios. We 
obtained correlated error processes by considering probe field propagation losses in a 
coherent feedback network, but analogous difficulties would result from any syndrome 
extraction mechanism in which the ancillary qubits used to accumulate the stabilizer 
generator values are subject to decoherence. For example, an idealized continuous-time 
measurement of the parity of a pair of register qubits Qi, Q2 in our framework [151 [6] can 
be thought of as corresponding to the standard quantum computational circuit diagram 
on the left, below: 



A 
Qi 
Q2 



H 



H 



A 
Qi 
Q2 



H 



E 



H 



The propagation losses we have considered essentially correspond to an error process 
acting on the ancillary qubit A, which takes place between the controlled-Z gates as 
indicated by the E gate in the above-right diagram. Of course, ancilla decoherence in 
quantum error correction has been treated in detail in the literature on fault-tolerant 
computing, as described for example in (TTJ HE] ■ In our context it is natural to assume 
that the dominant type of error process acting on the syndrome probe fields is optical 
loss, and in what follows we will show that it is possible to improve the robustness of 
the type of autonomous quantum memory we consider simply by optimizing the circuit 
layout. 

There turns out to be an interesting connection between circuit layout and 
robustness to propagation losses in our approach, for autonomous quantum memories 
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based on subsystem codes. Also known as operator quantum error correcting codes, 
subsystem codes are generalizations of decoherence free subspaces, noiseless subsystems, 
and quantum error correcting codes, which have gained popularity in recent years 
because of the large class of encoded logical operators these codes induce, which allows 
for simplified error recovery fZ3\ 121] . In [7] we considered a nanophotonic circuit for 
implementing a continuous-time version of the Bacon-Shor code, and noted that it is 
possible to reduce the circuit complexity by taking advantage of the subsystem flexibility 
in choosing register-qubit operations for corrective feedback. Here we further note 
that the subsystem structure also presents key advantages for syndrome extraction, as 
(following a fundamental insight discussed by Aliferis and Cross [12]) we can route the 
probe fields along paths that push the extra errors induced by optical propagation losses 
onto the unimportant gauge qubit degrees of freedom. Consequently, the correlated 
errors described above for the five-qubit code and in [8] for the nine-qubit code are no 
longer present and the storage fidelity of the encoded qubit is substantially increased. 
We will quantify the performance improvement using numerical simulations, below. 
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Figure 1. Left: Standard probe network for Bacon-Shor nine-qubit code. Right: 
Subsystem routing for optimal loss protection. 



To illustrate the robust routing strategy, we consider the Bacon-Shor six-body 
stabilizer generator Z^Z-jZ^Z^Z^Z^ which is one of two such Z-string operators for the 
bit-flip subnetwork El H5] . In our nanophotonic circuit model, optical propagation 
losses experienced by the probe field used to monitor this stabilizer generator will give 
rise to the following coupling terms, four of which are network-induced errors that the 
code is not designed to correct (while L n is a correctable error and L n+5 merely dephases 
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the stabilizer generator eigenstates, which does not compromise the quantum memory): 
L n = a6Z 8 , L n+ i = a6Z 5 Z 8 , 

L n +2 = o>6Z 2 Z 5 Z 8 , L n+ 3 = a6ZiZ 2 Z 5 Z 8 , 

L n +A = OidZ i ZiZ 2 Z^Z 8n L n+5 = a9Z 7 Z /L Z 1 Z 2 Z 5 Z 8 . (13) 

In deriving the above set of coupling terms we have assumed a geometrically simple 
routing of the probe beam, assuming the register qubit-resonator components are 
arranged in a 3 x 3 grid (which allows us to simplify the feedback signal routing as 
discussed in [7j), as shown in the left-hand schematic of Fig. [I] The scattering order is 
Z 8 — > Z 5 — > Z 2 — > Z\ — > Z 4 — > Z 7 , which bears a clear relation to the correlated error 
terms shown in Eqs. (13). However, we can measure the same stabilizer generator by 
routing the probe beam to the components in a different order as shown in the right- 
hand schematic of Fig. [TJ Z 8 — > Z 7 — > Z 4 — > Z 5 — >■ Z 2 — > Z 1 . In an idealized model with 
no propagation losses the scattering order should make no difference since the single- 
qubit Pauli operators commute. However, with propagation losses the second routing 
scheme leads to the following coupling terms in place of those of Eqs. (13): 

L n +i = <yOZ 8 , L n+2 = a8Z 7 Z 8 , 

L n +3 = olQZ^Z 7 Z 8 , L n+i = a8Z 5 Z 4i Z'i'Z 8 , 

L n +5 = a8Z 2 Z 5 Z 4 Z 7 Z 8 , L n+6 = a6Z 1 Z 2 Z b Z A Z 7 Z 8 . (14) 

In the subsystem structure of the Bacon-Shor code, the operators Z\Z 2 , Z 2 Z%, Z^Z 5) 
Z^Zs, Z 7 Z 8 , and Z 8 Zg are all logical operators on the unimportant gauge-qubit degrees 
of freedom. Consequently, each of the coupling terms in Eqs. ( 14 ) corresponds to either 
a single logical qubit error or the product of a logical qubit error and one or more 
gauge qubit errors. Since the single qubit errors taking place on the logical space are 
protected by the network, this implementation of the probe mechanism is substantially 
more tolerant of propagation losses, as we illustrate using numerical simulations in Fig. [2} 
In Fig. |2]we display the average fidelity decay F(t) = (^/ \p t \^ } over 10 4 simulated 
quantum trajectories (each) for lossy quantum memories implementing the five-, six-, 
seven- and nine-qubit codes J2H1 1221 [23] , which may be compared directly with analogous 
results from our prior work on other codes [HI [3, E]- In Fig. [3] we display the average 
decay (again over 10 4 quantum trajectory simulations for each code) of an alternative 
performance measure, 

F*(t) = maXi* e[tjt+r] F(t*). 

The quantity F* represents an easily computable, convenient statistic which we believe 
is a more meaningful measure of performance of a realistic quantum memory than the 
canonical fidelity measure. Our definition of F* is motivated by the observation that, 
in any realistic quantum memory, there must be a finite latency of error correction. 
This behavior is clearly illustrated for our class of models by the individual quantum 
trajectory simulations of fidelity versus time shown in Fig. |4} 

Because of the (random) time delay required for an error-correcting controller 
with finite-strength measurement and feedback to restore the register state after a 
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Figure 2. Decay of fidelity (^olptl^o) for the five-, six-, seven- and nine-qubit 
quantum error correcting codes with loss parameters 9 = {0, 1, 2.5, 5, 7.5, 10}7r/1000 
(top to bottom curves). For consistency with [5] the feedback strength 17 = is 
set to a constant value of 200 in each case. Note the substantial difference in the scale 
of the Y-axis for the bottom right corner plot — this is for the nine-qubit Bacon-Shor 
code with the optimally permuted probe network. 
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Figure 3. Averaged finite time-horizon fidelity for 10 4 quantum trajectories (each) 
of lossy five-, six-, seven-, and nine-qubit codes with loss parameter 9 — and 
t = 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35 and 0.4 (bottom to top traces in each sub-panel). 



decoherence event, fidelity does not decay monotonically along individual quantum 
trajectories [25]. This behavior has a pronounced effect on the appearance of a simple 
plot of average F(t) at small t as some trajectories in the ensemble will experience 
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Figure 4. Fidelity trace for single-shot quantum trajectory simulations of five-, six- 
seven-, and nine-qubit codes. 



errors at very early times without recovering immediately; this sub-ensemble induces 
the steep initial decay transient seen for example for t < 0.05 in the nine-qubit panel 
of Fig. [5J If we recognize that many (even most) of these trajectories will in fact regain 
F ~ 1 after a finite delay, as seen in the examples of Fig. [4], it seems intuitive to adopt 
a performance measure such as F*(t) that looks ahead over a window of time in each 
trajectory to check for such a recovery. Of course if an additional error should occur 
within a given trajectory before the feedback network has had time to correct the initial 
one, the encoded information can be lost and F*(t) — > accordingly. 
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